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THE UNUSUAL X-RAY MORPHOLOGY OF NGC 4636 REVEALED BY DEEP CHANDRA OBSERVATIONS: 
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ABSTRACT 

We present Chandra ACIS-I and ACIS-S observations (~200 ks in total) of the X-ray luminous 
elliptical galaxy NGC 4636, located in the outskirts of the Virgo cluster. A soft band (0.5-2 keV) 
image shows the presence of a bright core in the center surrounded by an extended X-ray corona and 
two pronounced quasi-symmetric, 8 kpc long, arm-like features. Each of this features defines the rim of 
an ellipsoidal bubble. An additional bubble-like feature, whose northern rim is located ~ 2 kpc south 
of the north-eastern arm, is detected as well. We present surface brightness and temperature profiles 
across the rims of the bubbles, showing that their edges are sharp and characterized by temperature 
jumps of about 20-25%. Through a comparison of the observed profiles with theoretical shock models, 
we demonstrate that a scenario where the bubbles were produced by shocks, probably driven by energy 
deposited off-center by jets, is the most viable explanation to the X-ray morphology observed in the 
central part of NGC 4636. As a confirmation to this scenario, radio jets extending towards the bubbles 
and a central weak X-ray and radio source are detected and are most likely the signs of AGN activity 
which was more intense in the past. A bright dense core of ~1 kpc radius is observed at the center 
of NGC 4636. A sharp decline in surface brightness from the core to the ambient gas is observed 
and is not accompanied by a variation in the temperature and thus could not be in thermal pressure 
equilibrium. However the bright core could be a long lived feature if the radio jets are acting as a 
balancing factor to thermal pressure or if the bright core is produced by steep abundance gradients. 
Subject headings: galaxies: individual (NGC 4636) - galaxies: ISM - galaxies: nuclei - X-rays: galaxies 



1. INTRODUCTION 

It is well known that active galactic nuclei (AGN) 
play an important role in the evolution of the hot 
gas in both individual galax i es an d clusters of galax- 
ies (e.g. IChurazqv et alJ 120001: iBriiggen fc Kaiser! 



2002! feavaliere et all I2002t iHoeft fc B ruggenl 1200 



Forman et alj|2005l IMcNamara fc Nulsedl200nV The so- 



called 'AGN feedback' can cause re-heating of the central 
cooling regions of a cluster and balance the cooling due 
to the X-ray emission (see IMcNamara fc Nulsenl 120071 
for a review) . AGN are also effective in shaping the mor- 
phology of the hot gas halos around individual galax- 
ies and clusters, e.g . giving rise to cavities and bubbles 
(|Birzan et al.ll2004l ). 

While several examples of subsonic bubble inflation due 
to AGN outbursts in individual gala xies are present in 
the literature (e.g . NG C 507, M 84; iKraft et al.ll200l 
iFinoguenov et al l 12008!) , supersonic bubble expansion, 
giving rise to shocks in the hot X-ray emitting gas, has 
only been o bserved in a handful of cases (e.g. M 87 , 
NGC 4552; iForman et all [20051 : iMachacek et all 120061) . 
Another possible case of supersonic bubble expansion 
could be that of NGC 4636, the dominant galaxy of 
a group on the outskirts of the Virgo Cluster (10° or 
2.6 Mpc on the sky to the s outh of M87, at a dis- 
tance to NGC 4636 of 15 Mpc; iTonrv et alJl200"l . This 



galaxy has been observed extensively by every major 
X-ray imaging observatory to date. Spectral imaging 
with Einstein showed the galaxy to be surrounded by 
an extensi ve halo of hot gas with a temperature of 
0.78 keV (|Forman et al.l Il985t ). Further studies using 
ROSAT and ASCA found gradients in both the tem- 
perature and abundance profiles , along with a very ex- 
tended halo (lAwaki et all 11994 iTrinchieri et al l 1199 



Matsushita et al.lll997t IFinoguenov fc Jonedl2000UBuot> 



200 



The metallicity is about 
(2') region, decreasing to 



20001 : lO'Sullivan et all 
solar in the central 9 kpc 
0.2- 0.3 solar at ~ 45 kpc (10') from the galaxy cen- 
ter (jMatsushita et alJll997h . More recently, NGC 4636 
has been observed by both Chandra and XMM-Newton. 
The initial Chandra ACIS-S observation showed an un- 
usual morphology in the core of the galaxy, most notably 
quasi-symmetric spiral arm-like features and a possible 
ring of high surface brightness emission, thought to be 
the pr oduct of shocks d riven by a previous AGN out- 
burst (|Jones et al.ll2002l ). As a confirmation to the out- 
burst driven shock scenario, a change of the gas tem- 
perature across one of the arm s was obser v ed as well. 
From the same Chandra data, lOhto et aT] (|2003h per- 
formed a detailed spectral analysis, suggesting the pres- 
ence of a cavity in the X-ray halo on the west side of 
the core, with AGN activity being the most likely cause 
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also in this case. In this scenario, plasma from the AGN 
jets may be able to displace the gas of the X-ray halo, 
leaving an apparent void as currently observed in the 
X-ray mo rphology. Using d ata fr om XMM-Newton and 
Chandra, lO'Sullivan et al.l (|2005h mapped the tempera- 
ture and abundance structure of the inner halo of NGC 
4636, confirming the detection of a cavity to the west 
of the galaxy and revealing the presence of a plume of 
cool metal-rich gas extending beyond the cavity to the 
southwest. Both cavity and plume appear to be the 
product of past AGN activity. Moreover, their abun- 
dance determination shows a depression in the central 
few kpc (~ 0.2 so l ar) in disagreement with the results of 
iMatsushita et alj (|1997t) . 

In this paper we take advantage of the long Chandra 
ACIS-I observations (~ 150 ks) combined with the pre- 
viously observed ACIS-S exposure (~40 ks) to analyze in 
more detail the features revealed by the previous Chan- 
dra and XMM-Newton observations. In particular we are 
able to study in more detail the morphology and temper- 
ature structure in the observed cavities and look for gra- 
dients in temperature and abundance to test the shock 
scenario invoked for their formation (|Jpnes "eTaLll200l 
lOhto et al.ll200l lO'Sullivan et~alll2005h . Moreover we 
will investigate with unprecedented detail the morphol- 
ogy and the properties of the hot gas in the galaxy core. 

At the distance of NGC 4636 (15 Mpc; iTonrv et all 
120011) . 1" w 70 pc. Throughout the paper, all the 
uncertainties are quoted at la (68%) for one inter- 
esting parameter, unless otherwise stated. The abun- 
d ance estimates are relative to the cosmic values given 
in lAnders fc Grevessd (|1989D . 

2. CHANDRA DATA PREPARATION AND ANALYSIS 

NGC 4636 was observed twice by ACIS-S and twice 
by ACIS-I. In our analysis, we use the longer ACIS-S 
observation (performed on 2000, January 26; ObsID: 
323) and both the ACIS-I observations, performed 
one immediately after the other on 2003, February 
14 (ObsID: 3926) and on 2003, February 15 (ObsiD: 
4415). The Chandra data analysis has been performed 
using CI AO v3.4, which applies the newest ACIS gain 
maps, time-dependent ACIS gain corrections, and 
ACIS charge transfer inefficiency corrections. The 
background light curve during the observation has been 
examined to detec t perio ds of high background following 
iMarkevitch et all (|2003h . A flaring episod (lasting 
~ 15 ks) is observed during the ACIS-S observation, 
reducing the usable exposure time to ~ 37 ks. No 
significant flaring episodes are detected during the two 
ACIS-I observations, giving a total good integration 
time for the ACIS-I detector of ~ 149 ks. 



2.1. Background Subtraction 

We use the "blank-field" observations, processed iden- 
tically to the galaxy observations (i.e. as described 
above) and reprojected onto the sky using the aspect 
information from the galaxy pointings. For the ACIS-I 
observations, the "blank-fields" are renormalized to the 
background of the observation using the ACIS-S2 chip. 
For the ACIS-S observation, the chip used for the renor- 
malization is ACIS-S 1 owing to the extended source emis- 
sion on S3. The energy band from 9.5 to 12 keV (mostly 



dominated by charged particles) is used to perform the 
normalization. 

We follow the procedure of IVikhlinin et all (|2005h . to 
check if the diffuse soft X-ray background could be an 
important background component in our observation. A 
spectrum from the source-free region of the S2 chip was 
extracted, the renormalized "blank-field" background 
was subtracted and the background subtracted spectrum 
was fitted in XSPEC vll.3.2p (in the 0.4-1 keV band) 
with an unabsorbed mekal model, whose normalization 
is allowed to be negative. The residuals are found to 
be negligible and consistent with zero in all three obser- 
vations, suggesting that no additional soft background 
correction is necessary for our datasets. 

3. RESULTS 

3.1. X-ray morphology and surface brightness 

Fig. [TJi shows the merged 0.3-2 keV Chandra image 
of all three observations (two ACIS-I and one ACIS-S 
observation). The image was "flat-fielded" using an ap- 
propriate exposure map weighted for the emission of a 
0.6 keV thermal plasma. The galaxy presents a very 
bright central core of radius ~ 1 kpc (15") centered on 
the nucleus. A lower surface brightness region surrounds 
the nucleus extending out to ~ 6 kpc (80"). Two quasi- 
symmetric arm-like features are embedded in this lower 
surface brightness emission. These f eatures were previ- 
ously reported by iJones et al.l (|2002t ) who analyzed the 
shorter ACIS-S observation. However, combining these 
data with the two ACIS-I observations we are able to 
observe that the south-western arm is clearly part of an 
X-ray cavity extending as far as ~ 9 kpc (2') from the 
center. The hint of a similar structure can be observed in 
coincidence with the north-eastern arm-like feature. It is 
not clear from the X-ray image whether there is a cavity 
in this case. However the cavities are more evident if we 
remove the contribution from the general diffuse emis- 
sion of the galaxy. Fig. [TJd shows the merged Chandra 
image in the 0.3-2 keV band where a /3-model fitted to the 
galaxy diffuse emission was subtracted. This processing 
allowed to highlight fainter structures. Indeed, both the 
south-western and the north-eastern cavities are clearly 
visible in this Figure. Moreover two additional quasi- 
symmetric cavities can be seen, just to the East of the 
nucleus and to the North- West of it, with the rims of the 
former being better defined than the latter. 
The radio data overlaid on the X-ray image (Fig. [2]) 
shows a correlation with the cavities. However, the ra- 
dio lobes do not fill the SW cavity, or break through the 
cavity boun daries as in other galaxi es presenting cavities 
(e.g. M 84. iFinoguenov et al.ll2008l ). However, the lobes 
detected at 610 MHz (Giacintucci et al. in prep.) with 
the Giant Metrewave Radio Telescope (GMRT) extend 
further than those observed at 1.4 GHz and also the NE 
cavity seems to be filled by plasma detected at 235 MHz 
with the GMRT (Giacintucci et al. in prep.). The ra- 
dio lobes are weak, having a luminosity at 1.4 GHz of 
~ 1.4 x 10 38 erg s _1 . However the plasma creating the 
cavities remains partially undetected by the radio obser- 
vations of NGC 4636 performed so far. 

Figure [3] shows Chandra images of NGC 4636 in two 
different bands: 0.3-0.9 keV and 0.9-1.3 keV. A 2" width 
gaussian smoothing was applied to them to enhance the 
features visible in the two images. The morphology is 
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quite different between them, with the X-ray arm-like 
features becoming prominent only at energies larger then 
0.9 keV, a possible indication that the temperature in 
these features is higher than in the surrounding hot in- 
terstellar medium of the galaxy. 

3.2. Global X-ray properties 

The surface brightness profile of NGC 4636 is shown 
in Figure [4j The emission does not show a central peak 
and is quite flat at the center, then it starts to decline 
monotonically at ~ 5" (0.35 kpc) from the center, show- 
ing another 'plateau' at 25" < r < 60" (1.75 kpc< r < 
4.2 kpc). Both these features can be noticed looking 
at the Chandra image (Fig. QJi). To study the radial 
spectral properties, we subdivided the emission from the 
galaxy into 25 annuli centered on the X-ray peak. We 
extracted a spectrum from each annulus (excluding the 
point sources) using specextract, which generates source 
and background spectra building appropriate RMFs and 
ARFs. The background is taken from the re-normalized 
blank field observations using the same region as the 
source. The spectra were analyzed using XSPEC vl2 
(|Arna ud 1996) and fitted w ith a single-temperature APEC 
model ([Smith et al.l[200lh . combined with the XSPEC 
pro j ct deprojection model. The free parameters in the 
model are the gas temperature kT, the gas metallicity 
Z and the normalization. The spectral band considered 
is 0.5-4 kcV. All the spectra were rebinned to have at 
least 20 counts per bin. We have measured Nh from the 
X-ray data, finding it consistent (within la) with the 
Galactic value along the lin e of sight, as der ived from 
radio data (1.9 x 10 20 cm' 2 ; iStark et al.lll992D . The de- 
projected temperature profiles of NGC 4636 are shown 
in Figure \5jpt,. No abundance profile trend was detected, 
we therefore decided to fix the abundance at Z = 0.5 
times solar. Th i s is in agreement with the findings of 
10' Sullivan et al.l (|2005l ). who also found no trend in Z(r) 
for the azimuthally-averaged data. Interestingly, they 
did however find a strong azimuthal abundance varia- 
tion. The temperature shows clear evidence for a decline 
in the center where kT ~ 0.5 keV, and increases toward 
the outskirts of the galaxy (kT > 0.8 keV). Several sub- 
structures, due to the complex X-ray morphology of the 
galaxy, are observed. The electron density profile de- 
rived from the deprojection (Fig. [5Jd) shows instead a 
strong central peak and a 'plateau' coincident with the 
one observed in the surface brightness profile located at 
2 kpc< r <4.5 kpc from the center. 

3.3. The X-ray Bubbles 

The complex X-ray morphology of NGC 4636 is mainly 
constituted by bubbles and cavities, whose properties 
cannot be determined with a simple radial analysis. 
Therefore, we studied the features observed in the X- 
ray image in detail, performing a spatially resolved spec- 
tral analysis on these regions, looking for evidence of the 
feedback from the central AGN in the surrounding hot 
interstellar and intergalactic medium. 

3.3.1. The SW Bubble 

The most prominent features observed in the Chandra 
image are clearly the two quasi symmetric 8 kpc long X- 
ray arm-like structures, already observed by iJones et al.l 



(2002) and by lO'Sullivan et~aTl (l200l . The SW arm 
(SW1) is however part of a cavity which extends at least 
~ 5 kpc in radius to the North (Fig. 0, although the 
northern boundary of the bubble (SW2) is not as sharp 
and well defined as the southern one. We performed 
a spectral analysis in a strip perpendicular to the SW 
arm-like feature, dividing the strip into rectangles ~ 0.5' 
wide (Fig. [7J, to look for variations in temperature or 
metal abundance. Although the metallicity does not 
vary significantly across the bubble rim, a sharp varia- 
tion in the temperature was clearly detected coincident 
with SW1 (Fig. [H]) showing a temperature decrease 
from kT ~ 0.75 keV to kT ~ 0.64 keV, well above 
the measurement errors (~ 0.01 keV). A temperature 
jump was not observed across SW2 most likely because 
of the complicated geometry of the X-ray emission. 
Indeed, SW2 seems to be partially embedded in another 
bubble-like feature located just North of it. 



3.3.2. The NE and the E Bubble 

The symmetry of the X-ray arm-like features is highly 
suggestive of the presence of a symmetric bubble NE of 
the nucleus of NGC 4636. However the southern rim of 
the bubble is not clearly visible and it looks instead to 
be embedded in another round shaped bubble located 
to the East of the nucleus. If we examine the surface 
brightness profile of the NE cavity we find a shape which 
is very similar to the SW bubble (Fig. [5^). Perform- 
ing a spectral analysis across the northern cavity rim 
NE1 (Fig. [7]), we also observe a temperature jump across 
the rim (Fig. HJa). However, the scenario in this part 
of the galaxy is more complex because of the presence 
of an additional feature just East of the nucleus with 
the shape of another cavity. This cavity looks less elon- 
gated than the other two cavities observed. Therefore, 
we preferred to analyze the temperature variation across 
its southern rim El (the northern boundary E2 could 
be contaminated by the presence of the NE cavity) us- 
ing partial annuli instead of rectangular regions (Fig. [7]). 
The surface brightness profile has a shape similar to the 
one across the NE cavity, while the temperature profile 
shows a temperature jump similar to the one observed in 
the other two cavities (Fig. [TO]). 

3.3.3. The origin of the cavities: a simple shock model 

The most likely scenario for the origin of the cavities 
observed in the X-ray morphology of NGC 4636 is that 
they were the result of successive outbursts of the cen- 
tral AGN. In this scenario the jet propagated rapidly 
from the center, creating a long thin cavity which then 
inflated in all directions. Perpendicular to the axis of the 
jet, the expansion has an approximate cylindrical sym- 
metry, motivating our choice to use rectangular regions 
and therefore assume cylindrical symmetry. The expand- 
ing cavities drive shocks into the surrounding gas. 

A 1-dimensional, cylindrically symmetric, time- 
dependent hydrodynamic model was used to investigate 
the properties of these shocks for the SW bubble. In this 
model the sound speed in the rclativistic gas that fills the 
cavity (the piston that drives the shock) is assumed very 
high, keeping the pressure in the piston nearly uniform. 
The ratio of the pre-shock pressure to the post-shock 
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pressure determines the strength of the shock. As a re- 
sult, the shock is weakest (slowest) in the region close 
to the AGN and fastest in the region farthest from the 
AGN. 

The undisturbed gas was assumed to be isothermal, 
with electron density n Cigas oc r~ n , where r is the dis- 
tance to the AGN. Corresponding to this, the surface 
brightness profile of the undisturbed gas has the form, 
SB oc p 1-2 '', where p is the projected distance from the 
AGN. The power law index r\ was determined by fitting 
the SB profile in the region adjacent to the SW arm, out- 
side the shock, giving r\ = 1.3. Taking a cut perpendic- 
ular to the radius at a distance, d, from the AGN, gives 
the electron density profile n c gas oc (d 2 + ra 2 )~' ) / 2 , where 
zu is the cylindrical radius measured from the axis of the 
lobe. Using the distance, d, corresponding to 1.2', this 
was the form adopt for the density profile of the undis- 
turbed gas in the cylindrical flow model. The initial gas 
distribution was made hydrostatic by imposing a static 
gravitational field on the model. Shocks were initiated 
by depositing a large amount of energy on the symmetry 
axis and their evolution was followed using the cylin- 
drically symmetric hydrodynamic code until the shock 
had expanded to the observed size. Surface brightness 
profiles for the model shocks were determined by pro- 
jecting the axially symmetric models onto the sky, using 
the density and temperature dependent 0.3 - 1.8 keV 
ACIS-I response. The temperature scale was determined 
by assuming a preshock temperature of 0.65 keV. The 
left panel of Fig. [8] shows the surface brightness profile of 
the SW shock, perpendicular to the radius, in the region 
shown in Fig. together with model surface brightness 
profiles for Mach numbers of 1.42, 1.72 and 2.04 (com- 
puted assuming 7 = 5/3 in the gas). Based on its surface 
brightness profile in the vicinity of the shocks, the Mach 
1.72 shock model provides the best fit. 

The temperature jump expected for a Mach 1.72 shock 
is a factor of 1.74. For a preshock temperature of 0.65 
keV, that would make the postshock temperature ~ 1.13 
keV. However, as observed, gas behind the cylindrical 
shock is projected onto the surrounding unshocked gas. 
The shocked gas also expands adiabatically behind the 
shock, reducing its average temperature. As a result, the 
observed temperature profile shows a considerably lower 
temperature jump than expected for a Mach 1.72 shock. 
Complicating matters further, fitted temperatures for 
multiphase gas are not read ily determined fr om the emis- 
sion measure distribution (Vikhlinin 2006). In view of 
this, projected emission measure distributions were de- 
termined for the Mach 1.72 model, assuming a preshock 
temperature of 0.65 keV. These were used to construct 
multiphase spectral models in XSPEC and the resulting 
models were fitted to a single temperature model to de- 
termine the observed temperature profile corresponding 
to the Mach 1.72 shock model. The results are plotted 
as stars together with the observed temperature profile 
in the right panel of Fig. [51 Within the errors, both the 
temperature and density profiles of this region are con- 
sistent with the Mach ~ 1.7 shock model. 

The calculations of the physical parameters derived 
from the model were performed for the SW bubble. The 
three cavities present similar temperature jumps and sur- 
face brightness profiles, so similar physical parameters 



are expected. From the hydrodynamic model, the age of 
the shock, i~2x 10 6 yrs, is the time it takes to expand 
to its observed size. The age is only sensitive to the shock 
strength and the temperature of the unshocked gas. This 
age is slightly shorter than the ratio of the shock radius to 
the present shock speed, because in the model the shock 
strength decreases with time (i.e. the shock expanded 
faster when it was younger). The total energy which 
produced the shock was 10 56 ergs, roughly equal to the 
enthalpy, H — 4pV — 7 x 10 55 ergs, calculated in the 
assumption that the bubble is predominantly relativistic 
(7 = 4/3). The average mechanical power required to 
produce the bubble equals to P mec h ~ 1-6 x 10 42 erg s -1 . 

The age of the cavity is shorter than what is usually 
observed in the so-called 'ghost' cavities where no ev- 
ident 1.4 GHz radio emission is observed to fill them. 
However, the cavities in NGC 4636 are not really 'ghost' 
cavities and a weak radio luminosity is ind eed expected 
from t he mechanical power derived above. iBirzan et al.l 
(2008) found indeed a correlation (although with a large 
scatter) between the cavity (jet) power Pj et and the ra- 
dio synchrotron power L radlo , of the form P jet ~ L^ adio . 
The cavities observed in NGC 4636 fall in the lower-lef t 
corner of the diagram in Figure 5 of IBirzan et al.l ([2008) , 
with a ratio Pjet/ ' ^ radio ~ 10 4 , similar to a few objects in 
the plot presenting a radio luminosity four order of mag- 
nitudes brighter than NGC 4636. There is a large scatter 
in the correlation, that is not well understood. Given the 
associated jet powers, the low level of radio emission from 
X-ray cavities in so me of these objects i s indeed surpris- 
ing. As discussed in IBirzan et alJ (|2008l ). the steep radio 
spectra in these systems cannot generally be explained 
by synchrotron aging or by adiabatic losses. Presumably 
some other loss mechanism is at work. This raises inter- 
esting questions on the n ature of the plasma filling these 
cavities. As stated e.g. in lJetha et al.1 (|2008f) . the nature 
of this material could be very hot, low-density plasma 
or could consist, at least in part, of a relativistic plasma 
containing relativistic electrons and magnetic field, sim- 
ilarly to that observed in the lobes of radio galaxies. In 
the case of NGC 4636, however, the complex X-ray mor- 
phology and the lack of low frequency radio maps (that 
will be presented in a forthcoming paper; Giacintucci et 
al. in prep.) does not allow at present to constrain the 
nature of the plasma filling the cavities (as it was done 
in e.g. IBirzan et aLll2008h iJetha et al.ll2008h . 

3.4. The nuclear 'point source 

The X-ray morphology of NGC 4636 is characterized 
by the presence of a dense core having a radius of ~ 1 kpc. 
A hard band image of the core (2-7 keV band) reveals 
the presence of several point sources. An X-ray point 
source is coincident with the radio nucleus (within 0.5") 
and may be X-ray emission from the supermassive black 
hole (SMBH) at the center of the galaxy. We extracted a 
spectrum of this point source fitting a power-law model 
to it in XSPEC. The spectrum has - 230 counts (ACIS-S 
plus ACIS-I) in total and the best fit model gives a power- 
law index T = 2.5to.t and a flux of 2.3±0.4 x 10~ 15 erg 
cm~ 2 s _1 and of 6.0±1.0 x 10~ 15 erg cm -2 s _1 in the 
2-10 keV and in the 0.5-10 keV bands, respectively. This 
corresponds to a total luminosity of 1.6±0.3 x 10 38 erg 
s _1 in the 0.5-10 keV band. This could be the current 
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luminosity of the AGN at the center of NGC 4636. How- 
ever, given the low luminosity, the X-ray emission could 
be also due to a low mass X-ray binary coincident with 
the radio nucleus. In the latter case, this value can be 
considered as a more stringent upper limit to the X-ray 
luminosity of the AGN at the center of the galaxy than 
the one set by lLoewenstein et al.l (|2001l ). 

To derive an estimate of the mass supply rate for ac- 
cretion on th e SMBH for NGC 4636, we can apply the 
iBondil (|1952f ) theory of steady, spherical and adiabatic 
accretion. This requires T and n at "infinity" , in practice 
near the accretion radius (that in the case of NGC 4636 is 
ta = GMbh/c 2 s ~ 8 pc, where c s is the speed of sound). 
Unfortunately the spatial resolution does not allow us to 
measure the temperature and the density at such a small 
radius. However, we can estimate the accretion rate, us- 
ing the temperature and the density profiles measured in 
in our radial analysis f§ 13. 2\ . The accretion rate is given 
by the formula 

M Bondl = 8.4 x 10 21 M 2 T 7 5 3/2 ni g a" 1 , 

where M s is the SMBH mass in u nits of 10 s Mp, (for 
NGC 4636 M BH ^ 7.9 x 10 7 M , IMerritt fe Ferraresd 
2001), To. 5 is the temperature in units of 0.5 keV and n% 
is the density in units of 1 cm~ 3 . Assuming that the den- 
sity n«r _1 , from the thermal component normalization 
of the spectrum in the inner 5" (~ 0.35 kpc), we obtain 
that the density at the accretion radius is n ~ 5.5 cm~ 3 . 
The temperature profile was instead fitted with a second 
order polynomial, giving a value at the accretion radius 
of kT ~ 0.49 keV. Using this values of n and kT, we 
find that M Bo ndi — 4.7 x 10 -4 M Q yr _1 . If at very small 
radius this gas joins an accretion disc with a standard ra- 
diative efficiency ij ~ 0.1, as in brighter AGNs, it should 
produce a luminosity: 

L acc = r]M Bondt c 2 ~ 2.7 x 10 42 erg s -1 , 

at least four orders of magnitude higher than the 
observed X-ray luminosity, pointing toward a highly 
inefficient accretion scenario. 



3.5. The bright core of NGC 4636 

The average density in the core is ^0.1 cm~ 3 and the 
thermal pressure (computed as 2.2nfcT, in the case that 
n = n p ) is ~ 9.7x 10 -11 dyn cm -2 . Several substructures 
are visible as well as can be seen from Figure [11] where an 
X-ray image of the central part of the galaxy is shown. 
This X-ray image (0.5-2 keV band) shows a cavity in the 
center of the core, although it is not clear, however, if we 
are dealing with a cavity or with a U-shaped enhance- 
ment in the surface brightness of the core. Interestingly, 
the position of this cavity is coincident with the inner 
radio jet detected at 1.4 GHz. 

We performed a finer spatial and spectral analysis of the 
core region, considering the X-ray emission west of the 
galaxy center (315°-405°), where a sharper decline in 
surface brightness is detected. This sector was chosen 
also because the discontinuity matches the curvature of 
a wedge centered on the nucleus. The surface brightness 
profile in this sector is shown in Figure [T2l In this profile, 
four regions with different behaviours can be identified 



easily. At the center a depression in the surface bright- 
ness (visible in Figure [TT]) is apparent. Outside the cen- 
ter (2" < r < 8") we observe a flat region, followed by a 
sharp decline at 8" < r < 18". At r > 18" the surface 
brightness flattens again then declines more gradually. 
We performed a spectral analysis in these four distinct 
regions to investigate variations in the temperature or 
in the abundance of the X-ray emitting gas. Because of 
the complexity of the emitting region and the number 
of counts detected we could not perform a meaningful 
deprojection, and we analyzed only the properties of the 
hot gas projected along the line of sight. A simple ther- 
mal model, plus a power-law to take into account the 
contribution of unresolved point sources, was considered 
(XSPEC model: wabs(apec+powerlaw)). The abundance 
was left free to vary in all the regions but the central 
(r < 2"), where it was fixed to 0.5 times solar because 
of the lower number of counts. Although a one tempera- 
ture fit does not always account for the complexity of the 
spectra (x 2 > 1.3), it gives an indication of the general 
temperature and metal abundance trend. Temperature 
and abundance profiles for the sector are plotted in Fig- 
ure [131 The temperature is constant at 2" < r < 18" 
(kT = 0.56 ± 0.01 keV) and it is slightly increasing right 
outside the core (kT = 0.60 ± 0.01 keV). The sharp de- 
cline in surface brightness (a factor of at least 5 over a 
range of 10") is indicative of a drop in the density of a 
factor of at least 2. The drop in density is clearly not 
accompanied by a corresponding variation in the temper- 
ature, thus if the thermal pressure is the only acting force 
in the central region of NGC 4636, the bright core is not 
in pressure equilibrium with the ambient gas and cannot 
be a long lived structure. At least two possible scenarios 
can be invoked to explain the presence of the bright core. 
One scenario involves the presence of the radio lobes 
which created the cavities observed in NGC 4636, where 
the pressure induced by the relativistic plasma could re- 
store the pressure equilibrium between the core and the 
ambient gas. In the other scenario the enhancement in 
brightness of the core region could be expl ained by an 
abun dance gradient (as in e.g. NGC 507; Kraft ct al. 
2004) . In Figure [131 a gradient in abundance between 
the flat region just outside the nucleus ( 2" < r < 8"), 
where Z — 1.0+J' 3 Z Q , and the adjacent regions, where 
Z ~0.6-0.7 Zq, is observed. This difference would not 
be enough to balance the pressure since the density is 
proportional to Z -1 / 2 , however the error bars are quite 
large and the fits demand for a more complex spectral 
model (1.3 < xt < 1-7). We performed a simple cal- 
culation to measure the diffusion timescale for Iron (and 
consequently for the other elements) to check whether the 
abundance gradient necessary to balance the pressure in 
the core could be long lived. For a thickness of the tran- 
sition region (between the peak in surface brightness and 
the flattening) of ~ 1 kpc, assuming Spitzer type colli- 
sion cross-sections for Fe, a density of 0.1 cm -3 , and a 
gas temperature of kT ~ 0.5 keV, the timescale for Fe 
diffusion is a few 10 13 yrs (i.e. several orders of magni- 
tudes longer than a Hubble time). For the lower atomic 
weight elements the diffusion timescale would be shorter, 
going approximately as the atomic weight cubed (e.g. a 
factor of ~ 40 times shorter for O than Fe). However, 
also considering lower atomic weight elements and all the 
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possible uncertainties in the measure of the abundance 
(mostly due to the strong correlation between kT, Z and 
flux in cool systems as NGC 4636), of the density and 
of the thickness of the transition region, the bright core 
could be very long lived if produced by an abundance 
gradient (and diffusion is covered by Coulomb collisions). 

4. CONCLUSIONS 

In this paper, we presented Chandra ACIS-I (~ 150 ks 
integration time) and ACIS-S (~ 40 ks) observations of 
the X-ray luminous elliptical galaxy NGC 4636, located 
in the outskirts of the Virgo cluster. The main results 
can be summarized as follows: 

• A soft band (0.5-2 keV) image shows the pres- 
ence of a bright core in the center surrounded by 
an extended X-ray halo and two prominent quasi- 
symmetric, 8 kpc long, arm-like features, both of 
which are rims of ellipsoidal bubble-like structures 
(NE and SW bubbles). An additional bubble- like 
feature (E bubble) is detected just south of the NE 
bubble. 

• We found that the temperature profiles across the 
bubble edges are sharp and characterized by the 
presence of temperature jumps of about 20-25%. 
On the other hand, the metallicity across the rims 
is constant. 

• Analyzing the core region of the galaxy we discov- 
ered the presence of a cavity around a weak two- 
sided radio jet (L ~ 1.4 x 10 38 erg s _1 ). Several 
point sources are detected in the hard band (2-7 
keV) in the core region. One of the point sources 
is coincident with the central radio source; its low 
X-ray luminosity of 1.6±0.3 x 10 38 erg s _1 may be 
consistent with that of a low mass X-ray binary 
coincident with a radio nucleus, however this lu- 
minosity may at least be considered as an upper 
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limit to the present X-ray luminosity of the central 
AGN. 

• We compared the observed temperature and sur- 
face brightness profiles across the bubble rims with 
numerical hydrodynamic shock models, demon- 
strating that a scenario where the bubbles were 
produced by shocks can explain the X-ray mor- 
phology. The shocks were probably driven by en- 
ergy deposited off-center by jets originating from 
the central AGN. A further confirmation of this 
scenario is the presence of radio jets extending to- 
wards the bubbles and of a central weak X-ray and 
radio nucleus, most likely the currently observable 
signature of AGN activity, which was more intense 
in the past. 

• For the SW bubble the shock has a Mach number 
of ~ 1.7, an age of ~ 2 x 10 6 yrs, and a total energy 
required to produce the shock of ~ 10 56 ergs. 

• A very bright dense (n ~ 0.1 cm~ 3 ) core of ~1 kpc 
(15") radius is observed in the galaxy center. The 
sharp decline in surface brightness (and conse- 
quently in density) in the core is not accompanied 
by a variation in temperature, not balancing the 
thermal pressure in the core with the ambient gas. 
This feature could be however long lived if the rel- 
ativistic plasma injected beyond the core by the 
radio jets is balancing the thermal pressure or if 
the surface brightness enhancement is produced by 
an abundance gradient. 
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Fig. 1. — Left: Chandra ACIS-I+ACIS-S image of NGC 4636 in the 0.3-2 keV band. The galaxy shows the presence of a bright central 
core (r ~ 1 kpc = 15"), surrounded by a lower surface brightness region extending as far as ~ 6 kpc (80") from the center. Two pronounced 
quasi-symmetric (8 kpc long) arm-like features define the rims of ellipsoidal bubbles. Also an additional bubble-like feature, whose northern 
rim is located ~ 2 kpc south of the north-eastern arm, is visible from the Chandra image. Right: Chandra ACIS-I+ACIS-S image after 
the subtraction of a /3-model fitted to the general diffuse X-ray emission. The rims of the ellipsoidal bubbles are more clearly visible and 
another bubble almost symmetric to the bubble-like feature observed south of the north-eastern arm, could be present north-west of the 
nucleus. 
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Fig. 2.— Chandra ACIS-I+ACIS-S image (0.5-3 keV) of NGC 4636 with the 610 MHz GMRT radio contours superimposed (from 
Giacintucci et al. in prep.). The radio contours shows a correlation with the cavities, however they do not fill them completely, or break 
through the cavity boundaries. The resolution of the radio image is 5.8" X 4.3", p. a. 48°, with the contours spaced by a factor 2 starting 
from ±3<r = 0.15 mjy/beam. An unrelated radio point-source is visible on the left of the nucleus, just below the NE cavity. 





Fig. 3.— Left: 0.3-0.9 keV Chandra ACIS-I+ACIS-S image of NGC 4636 smoothed with a 2" width gaussian; Right: 0.9-1.3 keV Chandra 
ACIS-I+ACIS-S image of NGC 4636 smoothed with a 2" width gaussian. The X-ray arm-like structures become prominent only in the 
higher energy band. In the figure 1" corresponds to ~ 4.2 kpc. 
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Fig. 4. — Chandra ACIS-I+ACIS-S surface brightness profile of NGC 4636 extracted in the 0.5-3 keV band. A central surface brightness 
peak is not present, being quite flat at the center. The profile declines monotonically at 5" < r < 25" (0.35 kpc < r < f .75 kpc), showing 
another 'plateau' at 25" < r < 60" (f .75 kpc < r < 4.2 kpc). 
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Fig. 5. — Left: Deprojected temperature profile for the X-ray gas observed in NGC 4636. A decline of the temperature is observed at 
the center, where kT ~ 0.5 keV, while the temperature is increasing toward the outskirts of the galaxy (kT > 0.8 keV). Right: Deprojected 
electron density profile for NGC 4636. The density shows the presence of a central peak and of a 'plateau' coincident with the flattening 
observed in the surface brightness profile, located at 2 kpc< r <4.5 kpc from the center. 
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Fig. 6.— Chandra ACIS-I+ACIS-S image of NGC 4636 (0.5-3 keV band). The three bubble-like features detected are labelled and 
identified by cyan ellipses. The yellow arrows are pointing toward the detected rims of the bubbles. 
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Fig. 7. — Chandra ACIS-I+ACIS-S image of NGC 4636 in the 0.5-3 keV band. The regions used to extract the spectra for the analysis 
of the three bubble-like features are shown and labelled in white. 
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Fig. 8. — Left: Surface brightness profile perpendicular to the shock front for the SW bubble. The three colored solid lines represents 
the prediction from a numerical hydro-dynamical shock model at different Mach numbers for the shock. The best fit shock model to the 
observed data has M = 1.72. Right: Temperature profile across the southern rim of the SW bubble. The temperature jump is consistent 
with the predictions of a shock model with M = 1.72. The zero points in r are coincident in the two panels and the values of r increase 
toward the South. 
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Fig. 9. — Left: Surface brightness profile perpendicular to the shock front for the NE bubble. Right: Temperature profile across the 
northern rim of the NE bubble. The shape of the surface brightness profile is very similar to the one observed in the SW bubble, moreover 
a temperature jump is observed in coincidence with the northern cavity rim (NE1). The zero points in r arc coincident in the two panels 
and the values of r increase toward the South. 
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Fig. 10. — Left: Surface brightness profile across the southern rim for the East bubble. The dashed vertical lines delimitate the regions 
used for the spectral analysis. Right: Temperature profile across the southern rim of the E bubble. The shape of the surface brightness 
profile is analogous to the one observed in the other two bubbles and a temperature jump is observed in coincidence with the southern rim 
(El) at ~ 2.5 kpc from the center of the bubble. 
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Fig. 11. — Chandra ACIS-I+ACIS-S image of the core of NGC 4636 (0.5-3 keV band). A cavity in the galaxy core is well visible and it 
is surrounded by a U-shaped brighter rim. 
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Fig. 12. — Surface brightness profile in a sector west of the nucleus (315°-405°), extracted in the 0.5-3 keV band. A depression in the 
center is well visible, followed by a flat region and then a sharp decline of the surface brightness. The dashed vertical lines delimitate the 
regions used for the spectral analysis. 
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Fig. 13. — Left: Temperature profile in the 315°— 405° sector. Right: Abundance profile in the 315°-405° sector. Both the temperature 
and the abundance profile are quite flat and do not present evidence of strong gradients with the radius. 



